Electron corrected Lorentz forces in solids and molecules in magnetic field 
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We describe the effective Lorentz forces on the ions of a generic insulating system in an magnetic 
field, in the context of Born-Oppenheimer ab-initio molecular dynamics. The force on each ion 
includes an important contribution of electronic origin, which depends explicitly on the velocity of 
all other ions. It is formulated in terms of a Berry curvature, in a form directly suitable for future first 
principles classical dynamics simulations based e.g., on density functional methods. As a preliminary 
analytical demonstration we present the dynamics of an H2 molecule in a field of intermediate 
strength, approximately describing the electrons through Slater's variational wavefunction. 
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Understanding the behavior of matter in large mag- 
netic fields is important both conceptually and practi- 
cally. The effect of a magnetic field on a non-magnetic 
(molecular, fluid or solid) electronic insulator is twofold. 
The first is on the electronic states, via a field induced 
modification of electron quantization. The states of free 
electrons in a magnetic field are split into Landau levels; 
in a general insulator, the electronic states or bands will 
turn into some generalized form of Landau levels, whose 
splittings are important when at large fields they become 
comparable with band gaps. These electronic effects can 
be efficiently calculated with modern density functional 
(DFT) methods well described in the literature. [I]. The 
second type of field effects, important in dynamics, are 
those on nuclear motion, via Lorentz forces. In an insu- 
lating or close shell system, where electron motion can be 
decoupled adiabatically, the effective Lorentz force F n on 
any ion n still depends on the electronic states. Schemat- 
ically, the total effective Lorentz force may range between 
the bare Lorentz force F n = (Q n /c)V x B (if electrons 
could be ignored) to zero (if electrons could be assumed 
to be tightly bound to nuclei, thus neutralizing them). 
Even if the effect on the electronic structure is very small 
and probably negligible for solid state applications, the 
electronic contribution to the Lorentz force is not negli- 
gible at all, but is so far without a clear understanding 
of how one may, at least in principle, calculate it. 

A pioneering fully quantum-mechanical treatment of 
nuclei and electrons in a magnetic field was formulated 
long ago, establishing a clear basis of principle. Q, 0| ■ The 
separation of the center of mass nuclear motion gives rise, 
apart from a trivial motional Stark effect, to a number 
of mass-correction terms in the electronic and nuclear 
hamiltonians, which provide additional couplings of the 
motion of different nuclei. Unfortunately, this approach 
is computationally very demanding and can be applied 
at best to small molecules or clusters. There exists to 
our knowledge neither an explicit formulation that can 
be used right away for ab-initio simulations, for example 
of DFT type, nor a direct application to real systems, 



beyond the hydrogen atom 0]. 

In this letter we pursue a formulation of electronic 
Lorentz forces based on the Born-Oppcnhcimcr approxi- 
mation, in principle suitable for state-of-the-art first prin- 
ciples simulations. As an illustration, we apply it to the 
dynamics in a magnetic field of the H2 molecule described 
by a simple variational wavefunction. Following 0, 
we will use the electronic Berry phases as a tool to cal- 
culate electronic Lorentz forces acting on the ions in a 
general insulating system. The BO approximation, which 
amounts to assume fixed and infinite mass nuclei, is made 
before the separation of the center of mass motion. By 
these assumptions we neglect the motional Stark effect 
and the mass correction terms. The former is well jus- 
tified as long as the particle speed is small compared 
to the speed of light, as in most condensed systems Q. 
Mass correction terms become important for exceedingly 
strong magnetic fields, when the insulating character may 
be lost and/or particles may behave relativistically, in 
which case these corrections should of course be included 
in the nuclear and electronic hamiltonian. 

We assume the ground state electronic wavefunction to 
be given, at each time step, by an electronic structure cal- 
culation, such as Hartree-Fock 0], or density functional 
theory (DFT) methods [8|. The nuclear hamiltonian of 
a generic insulating system in the BO approximation is 



Hn = 



A{R n ) - Xn 



U({R}). 



(1) 

The first term is the ion kinetic energy (P n = — iVi^); 
the second is the ground state expectation value of the 
electronic hamiltonian, including the magnetic field effect 
on the electronic structure. Xn appearing in the first 
term of eq. ([I]) is the so called Geometric Vector Potential 
(GVP) or Berry connection, given by 0, 0| 



(2) 



where ip c \({R}) is, in full generality, the many-body elec- 
tronic wavefunction in presence of the field, normalized to 
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the number of electrons: (ip c \\tp c \) = N e \, and depending 
parametrically on all nuclear coordinate. Whenever this 
wavefunction can be chosen real and single-valued, the 
GVP of eq. @ vanishes. 0] In a magnetic field, the wave- 
function cannot be made real, and the GVP is nonzero. 
After integration of the electronic degrees of freedom, 
a gauge potential associated with the GVP appears in 
the nuclear hamiltonian [1, H, @| . This gauge field plays 
the role of an additional magnetic field, one that couples 
only to the kinematic degrees of freedom, and not e.g., to 
the nuclear moments, which instead experience the real 
field. We restrict here for simplicity to large gap insula- 
tors, where the adiabatic approximation is well justified, 
and one can safely ignore all excited electronic states in 
the expression of Xn- In that case moreover the spin 
susceptibility is minuscule, and spin effects can also be 
neglected. 

We derive from eq. |T]) the nuclear equations of mo- 
tion (EOMs) from the Heisenberg time evolution of the 
positions and velocity operators. Classical equations of 
motion are then obtained through Eherenfest's theorem 

M n V n = -VH^+tQJcjKxB+^KnxOW (3) 

TO 

Eq. {3]) resembles that of a charged ion in a field, but for 
the last term, which is precisely the electronic Lorentz 
force. The gauge invariant quantity J~2( nm ) is the Berry 
curvature which plays the role of an effective magnetic 
field (gauge field) in the 3N— dimensional space spanned 
by the ionic degrees of freedom (N is the number of nu- 
clei). It is given by 

n (nm) = _ 2 Im (V^cl | x | Vj^Vel) (4) 

and has the dimensions of a magnetic field, in fact propor- 
tional to the external field (B/c) when sufficiently weak. 
Unsurprisingly, the force on ion n now depends upon the 
velocity of all other ions m through the off-diagonal terms 
in eq. . If the electrons were infinitely tightly bound to 
the nuclei (or equivalently, in the limit of large separation 
between the ions), these off-diagonal terms would vanish. 
In that regime, ip c \ collapses to a sum of products of sin- 
gle particle orbitals centered around the nuclei. Each 
is rapidly decaying in space, and each is dragged rigidly 
along by its nucleus as it moves, thus providing total mag- 
netic screening for every ion. Reality is of course very far 
from that limit, magnetic screening is only partial, and 
must be calculated explicitly. 

The basic ingredient for computing the Berry cur- 
vature eq. (j4|) is the derivative of the electronic wave- 
function with respect to atom position R n . That can 
be obtained in an electronic structure calculation by fi- 
nite differencing the electronic wavefunction ip e \ for two 
atomic configurations, compensating the arbitrary phase 
of the wavefunctions (covariant derivative). Alterna- 
tively, |Vij n Vei) can be obtained by linear response to 
an atom displacement. 



In order to provide a first exemplification of the elec- 
tronic Lorentz forces, with a direct analytical and quan- 
titative insight into the properties of the Berry cur- 
vature eq. ([3]), we consider here as a simple example 
the classical dynamics a neutral homonuclcar diatomic 
molecule in a field. If we set the molecule in motion and 
freeze the vibrational and rotational degrees of freedom 
(Vi = V 2 = V), the EOMs are 

MV = (Q/c)V x B + V x (W 11 ' + ft (12) ) (5) 

(inversion symmetry requires fi^ 11 ) = f2( 22 ) and fi^ 12 ) = 
il*- 21 )). Since the molecule is neutral, the total Lorentz 
force must vanish and f^ 11 ) + Jl^ 12 ' = —Q(B/c) must 
screen completely the nuclear charge Q. In the het- 
eronuclear diatomic case, the nuclei are not screened in- 
dividually, but only globally, in the form ~J2 nm fl^ nm ^ = 
— (B/c) ~J2 n Q n [10]. If we consider a rigid rotation of the 
frozen molecule, then V\ = — V 2 = V around the center 
of mass and the EOM is 

MV = {Q/c)V x B + V x - fJ (12) ) (6) 

The electronic screening field fi r = fi^ 11 ) — fi^ 12 ) depends 
upon the internuclear distance and is related to the elec- 
tronic Berry phase 7 accumulated on a single period of 
rotation: 

7 = I x(R) -dR= II n T {R) ■ n dS (7) 

JdC(R) JJc{R) 

where C(R) is a circle of diameter R R being the inter- 
atomic distance. As shown in Ref. [ll|, 7 is the electronic 
contribution to the rotation- induced magnetic moment. 

We now wish to address the generically rotating, trans- 
lating and vibrating molecule. In general, that calcu- 
lation can be done by first principle implementation of 
eq. (J3j> and ((4]). In order to make the illustration more 
explicit while keeping it simple, we will focus on a hy- 
drogen molecule, where we can obtain essentially analyt- 
ical results by describing the electronic structure through 
Slater's variational LCAO approximation based on two 
Is orbitals [13]. To include the field, we form linear com- 
bination of gauge-including atomic orbitals (GIAO) 

tpei = ci (pi(r - i?i) + c 2 &2(r - R 2 ) 

<f n = exp[-(ie)/c$(Rn -> r)} ip ls (r - R n ) (8) 

where the phase factor $(r — > r') is the integral of the 
vector potential along the line connecting the points r 
and r' [131 ] . ipi s is the hydrogenic Is radial wavefunc- 
tion for a nuclear charge a, taken as variational parame- 
ter [13]: ipis(r) — (a 3 /^) 1 / 2 e~ ar . For a fixed value of the 
internuclear separation, a minimizes the sum of the elec- 
tronic energy plus the nuclear-nuclear coulomb repulsion, 
and varies from 1 in the limit of large separation to ~ 2 in 
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FIG. 1: Top: SDf' (solid line) and n[ 12) (dashed line) in 
units of (B/c) as a function of the interatomic distance, for 
the H2 molecule. Bottom: The effective magnetic field felt by 
each proton 1 + $"l r , in units of (B/c). Left panels, field nor- 
mal to the molecule axis; right panels field parallel molecular 
axis. The vertical line indicates the theoretical equilibrium 
distance. 



the limit of small separation where the H2 molecule col- 
lapses to a He atom. We neglect the electron spin due to 
the large gap between the singlet and triplet (~ 10 eV). 
Owing to its simplicity, this trial wavefunction gives for 
H2 a crude equilibrium distance of 1.00 A (experimental 
0.74 A) and a barely reasonable dissociation energy of 
4.235 eV (experimental 4.476 eV). However, it illustrates 
the electronic Lorentz forces very well. At weak field, the 
one electron matrix elements are 



H nm = H nm exp [-(ie/c) $(R m -> Rn)} 
S„ m = S nm exp [-(ie/c) $(R m -> Rn)] 



(9) 



where H nm = {ip n \ H° | tp m ) and S nm = {ip n \ <p m ) are 
the field-free matrix elements. 

We computed fl^ 11 ^ and JV 12 ' with the trial wavefunc- 
tion of eq. dU), directly from eq. . The resulting ex- 
pression, analytic, contains a number of terms and is too 
long to be shown here. Since GIAOs give the first order 
correction in (B/c) to the zero field wavefunctions, we re- 
tain only the linear (B/c) term in the Si's, where the zero 
order term is zero. These quantities depend on the inter- 
atomic distance and on the angle between the molecule 
axis and the magnetic field, preserving the cylindrical 
symmetry of the system, independent of the gauge. In 
the following J~2( nm ) will indicate the coefficient of the 
first order term in (B/c). 

The upper panels of Fig. Q] show jV 11 ) and as 
a function of the interatomic distance, for two relative 
orientations of the molecule with respect to the magnetic 
field, which we choose parallel to the z axis. For large in- 
teratomic separation R, i.e. in the dissociation limit, the 



off-diagonal JV 12 ) vanishes and f^ 11 ) reaches the asymp- 
totic value —1. This reflects the fact that the motion of 
the two protons is decoupled in the dissociation limit and 
the electronic screening of the individual nuclei is com- 
plete. For small interatomic separation, both fi^ 11 ) and 
fi (12) (and by symmetry n ( - 22 '> and J^ 21 )) tend to -1/2, 
recovering the correct screening of an isolated He atom. 
For arbitrary interatomic distance, the sum of fl( n > and 
fj( 12 ) is identically —1, which fulfills eq. §5§ and warrants 
total screening of the center of mass motion. The lower 
panels of Fig. Q] show the effective relative magnetic field 
at the proton site, (1 + fi r ) as a function of the inter- 
atomic distance, for two orientations of the molecule rel- 
ative to the field. At the equilibrium distance of R ~ 1 A, 
the proton's Lorentz force is that of a reduced effective 
charge ranging between +0.4|e| and +0.6|e|, depending 
on orientation. (The true reduction is actually weaker, 
since our variational wavefunction slightly overestimates 
the electronic screening). The approximate rotational g- 
factor [3] computed by eq. (J7J| in the present approxima- 
tion is 0.62, in fair agreement with experiment [14[, and 
with more accurate calculations yielding 0.88. This 
sort of error does not impair the value of the present ap- 
proximation as an analytical illustration of the method. 

Armed with the Berry curvatures J7(™ m ) - now known 
analytically - and with a simple paramctrization of the 
interatomic potential U(R), we can describe the classi- 
cal dynamics of H2 in a field. As in a first principles 
molecular dynamics simulation, we integrate the equa- 
tion of motion eq. j3j for the H2 molecule, exploring the 
effect of different initial (Cauchy) conditions [HI]. By 
way of example, we start by compressing or stretching 
the molecular bond and let it free to vibrate/rotate. The 
restoring force sets the two protons in motion initially in 
the radial direction but soon their trajectory is deflected 
by the Lorentz forces. If the molecule initially lies in the 
plane perpendicular to the magnetic field, the resulting 
orbits resemble cycloids in that plane. The resulting tra- 
jectories are shown in the upper panels of Fig. f2J The 
sense of rotation is determined by the initial condition of 
stretching or compression. 

To clarify the effect of the electronic Lorentz forces 
on the EOM we show in the lower panels of Fig. [5] the 
same trajectories now obtained by neglecting the Berry 
curvatures from eq. ([3]) - i.e. only retaining the bare 
Lorentz force (Q/c)V x B. When the Berry curvature 
is included, the angular velocity of the cycloid is reduced 
by the screening action of the electrons by a factor ~0.4, 
which is also the average fraction of magnetic field felt 
by the ions during the vibrations around the equilibrium 
position (see r2 r in Fig. [1]) . We note that this reduction 
factor measures the strength of the effective magnetic 
field at the proton site, and differs from the rotational 
g-factor (here 0.62), which measures instead its integral 
over the orbit spanned in a full rotation. Summing up, 
implementation of eq. ([3]) and ^ yields a description of 
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FIG. 2: Simulated trajectories of the H2, in presence of mag- 
netic field. Left panel: neglecting the Berry curvature term. 
Right panel: with the Berry curvature term. In both cases, 
the molecule was set initially parallel to the x axis and the 
bond compressed 0.1 A. The molecule rotates in the counter- 
clockwise sense. For a field 10 T, the estimated rotation pe- 
riod is of 2-10~~ 5 s. The units of the plots are A. 



nuclear motion whose accuracy is only limited by that of 
the underlying electronic calculation. 

All of the above is classical nuclear motion. In the 
quantum EOMs for the nuclei, the presence of the GVP of 
Eq. ^ must be considered. In particular, starting from 
the zero field nuclear hamiltonian of a diatomic molecule, 
we may explore how the GVP term of Eq. @ influences 
the roto-vibrational spectrum. In the symmetric gauge, 
A(x) = (1/2)1? x x one can separate center of mass 
and relative distance in Eq. (JTJ) . The distance hamil- 
tonian is H R = 1/(2//) [p R + (e/c)A(R) + xr? + U(R) 
where xr = (e/2c)B x R(S - 1)/(S + 1), S being the 
overlap integral between the two atomic orbitals. The 
curl of Xr is O r . For ordinary laboratory magnetic 
field intensities, the field can be considered as a per- 
turbation, therefore we expand the hamiltionian up to 
first order in the field H = Hq + H\. Hq is the un- 
perturbed hamiltonian of a harmonic vibrating rotator, 
Hi = (eB) / (2fj,c)L z 2S/ (S +1) is the rotational paramag- 
netic term. The basis of the unperturbed hamiltonian Ho 
can be labeled by the quantum numbers (n, I, m), and the 
spectrum is given by the rotovibrational levels of the di- 
atomic molecule. To first order in (B/c), Hi removes the 
degeneracy of the rotational levels is removed according 
to the usual Zeeman splitting E n \ m — E^ m +gRyL nuc B m, 
where gn is the rotational (/-factor. The eigenstates have 
corrections of order (B/c) due to the term S/(S + 1), 
which mixes an unperturbed state with other states with 
the same (l,m) and different n. However, for a field of 
1 T, the coupling between vibrational states is of the or- 
der of 10 -5 cm -1 , much smaller than the centrifugal and 
anharmonic couplings, which are of the order of 10 cm . 
To second order in (B/c), there will be diamagnetic shifts 
affecting mainly the states with L z = 0. Despite the rota- 
tion, the canonical angular momentum L z is a conserved 
quantity. In the symmetric gauge, the mechanical angu- 



lar momentum is 

where / and uj are respectively the momentum of iner- 
tia and the angular velocity. The second term in the 
right-hand side is always positive, depends on the vibra- 
tional state of the molecule, and is of the order of (B/c). 
The physical rotation of the molecule correspond to the 
expectation value of Iui; when m =/= the difference be- 
tween Iu> and L z is negligible because it is of order (B/c). 
However, eq. (fT0|) shows that even when L z — 0, a small 
amount of rotation still exists (see Fig. [2]). 

In conclusion, we presented a convenient formalism to 
calculate the all important electronic contribution to adi- 
abatic Lorentz forces for atomistic dynamics in a mag- 
netic field, based on Berry connections and ideally suit- 
able for future ab-initio simulations. We demonstrated 
its validity and applicability in the simple example of H2 
where variational wavefunctions provide approximate but 
analytical results for the the Berry connection and Berry 
curvatures. The example demonstrates a weak field in- 
duced coupling between rotations and vibrations. This 
work is supported by MIUR COFIN No. 2004023199- 
003, FIRB RBAU017S8R operated by INFM and MIUR 
FIRB RBAU017S8. 
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